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ABSTRACT 

Fuzzy logic has proven to be a simple and robust method for process control. Instead of requiring a complex model 
of the system, a user defined rule base is used to control the process. In this paper the principles of fuzzy logic control are 
applied to Particle Tracking Velocimetry (PTV). Two frames of digitally recorded, single exposure particle imagery are used 
as input The fuzzy processor uses the local particle displacement information to determine the correct particle tracks. Fuzzy 
PTV is an improvement over traditional PTV techniques which typically require a sequence (>2) of image frames for 
accurately tracking particles. The fuzzy processor executes in software on a PC without the use of specialized array or fuzzy 
logic processors. A pair of sample input images with roughly 300 particle images each results in more than 200 velocity 
vectors in under 8 seconds of processing time. 


2. INTRODUCTION 

The control of complex processes has been aided by the development of fuzzy control systems. Fuzzy control 
systems have been used to control traffic flow, appliances, and even subways for optimal energy efficiency and passenger 
comfort. 1 Fuzzy logic employs a tolerance for imprecision to achieve system control. 2 An exact model for the system inputs 
and outputs is not required. Fuzzy inference control utilizes membership functions and a rule base developed by the user to 
process information. The physical mechanisms underlying the process are impertinent to the controller. The process of 
identifying and tracking particles in a flow is a good candidate for fuzzy control since the procedure is not clear cut, but 
involves some gray area decisions. 

Particle Tracking Velocimetry (PTV) is a planar velocity measurement technique whereby the positions of particles 
entrained in a fluid are recorded over a planar region at more than one instance in time. 3 A light sheet is used to illuminate 
the particles in the fluid. The exposures can be recorded either on a single flame, or on independent frames. When 
independent frames are used, the velocity vector direction information is easily determined. However, mechanical shifting 
between exposures recorded on a single frame also reveals direction information. 4 

Double exposure, single frame imagery are more commonly processed via correlation techniques and require high 
data densities. 4 The correlation techniques rely on the local average particle displacement to determine the net velocity over 
the interrogated region. PTV techniques typically require more than two exposures and can perform well in both high and 
low data density cases. 5,6 The PTV techniques require more time history information in order to make an accurate 
identification of a valid particle track, which requires more data storage and/or higher framing rate cameras. 

This paper describes the application of fuzzy inference techniques to track individual particles recorded on a pair of 
single exposure frames. PTV is a good candidate for fuzzy control since the particle trajectories are never identical and the 
velocity field varies spatially across the image plane. Local particle displacement information is used to identify candidate 
particle tracks. The fuzzy inference processor is used to determine the most probable particle trajectories based on common 
sense rules that an observer would use to identify particle tracks. This is the first known application of fuzzy control to PTV. 
Similar work has been done in applying neural networks to particle tracking. 7 

A sample flow field is recorded and reduced using the Fuzzy PTV reduction technique. The fuzzy logic PTV runs in 
software on a PC, without specialized array processors. A sample data set was reduced in under 8 seconds yielding over 200 
velocity vectors. 
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3. PARTICLE TRACKING VELOCIMETRY 


The PTV technique described herein employs direct digital recording of the particle images onto a CCD array. Two 
successive single exposure video fields are acquired from the CCD camera. The two image fields are independently processed 
to determine the particle image centroids. The particle image amplitude and shape information are used to compute the 
intensity weighted mean particle location. The centroid locations are stored and used to determine candidate displacement 
vectors. Velocity vectors are calculated by dividing the displacement vector by the inter-exposure time. 

The particle centroids on the initial Frame #1 are used as starting points for possible particle displacements. The 
user specifies a search region radius, R„ typically 10-20 pixels, to search for Frame #2 particles. Each Frame #2 particle 
within a radius R, from the initial particle centroid is a candidate displacement vector. All possible displacements of the 
ini tial particle to the second particle locations within the search region are recorded and stored as lists of candidate 
displacement vectors for that initial particle. Hence, for high data density areas, many initial particles may be competing for 
the same second exposure particle centroids. At this stage in the processing, the vector field is very convoluted and noisy. 

The fuzzy inference processor operates on these lists of candidate displacement vectors to determine the most likely 
displacement vector for each initial particle centroid location. 

The list of vectors for each initial particle is compared to all other initial particle displacement vector lists to 
determine if there is any commonality. If two separate initial particles do claim the same second particle, then all possible 
vector pairs between these two initial particles are compared. The main assumption is that if two initial particles are close 
enough to interact (claim the same second particle) then the pair of vectors that look the most similar (in direction and 
magnitude) must be the correct pair of displacement vectors for the two separate initial particles. This assumption also holds 
for tertiary and higher interactions. 


4 FUZZY LOGIC PROCESSOR APPLIED TO PTV 

There are four inputs to the Fuzzy PTV processor for each vector pain distance between the vector midpoints in 
pixels (Sep); average vector magnitude (Mag); difference in vector magnitudes (MagDif); and the sum of the squares of the 
differences of the x- and y-components of the two velocity vectors (Delta). The first and third measures, when combined, act 
as a velocity gradient measure. The last measure operates in the opposite manner to the dot product. This measure is small 
when the vectors are similar in magnitude and direction and large when they are different. Each input measure is assigned to 
a fuzzy set, where the degree of membership for each element in the set varies between 0 and 1. Standard 25-50% 
overlapping triangular input membership functions are used as shown in figure 1. The degrees of membership for each input 
are processed through a rule base of "IF...THEN" blocks. The rule base defines an output fuzzy set For a given vector pair, 
up to 16 rules may fire depending on the number of unique combinations of membership values. In lieu of the more common 
centroiding technique, the Fuzzy PTV processor output is computed via the singleton technique with a weighted average, 
which is computationally simpler. 8 

The Fuzzy PTV rule base is shown in matrix format in table 1. Each input measure has three elements in its fuzzy 
set: {Small, Med, Large}. Every fuzzy set has a corresponding set of membership functions. The membership functions 
assign a numeric weight, or degree of membership |x, for all elements in each input measure's respective fuzzy set. The 
number of non-zero memberships depends on the amount of overlap of the membership functions. The midpoints of the 
membership functions are: Sep_Mid=R J ; Mag_Mid=Rj/2; MagDif_Mid=5 pixels; Delta_Mid=5 pixels. R s is the search 
region size in pixels, which is nominally 20 pixels. 

The input measures are computed for each vector pair analyzed. The fuzzy processor output from the rule base is 
the fuzzy set of confidence in the vector pair. The set of confidence singletons and their associated numerical weights are: 
{High=0.8; Med=0.5; Low=0.1; None=0). For example, in table 1 with Sep=Small, Mag=Med, MagDif=Small, and 
Delta=Med yields Conf_Out=Med. 
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Table 1: Rule base for Fuzzy PTV processor. Each unique combination of inputs defines an output 
element of Conf^OuL 

All unique combinations of the degrees of membership for each fuzzy set for each input measure are combined via a 
fuzzy "AND", which is the minimum value of the four membership values. The actual confidence level in the vector pair is 
computed by taking the product of the fired rule output value (Conf_Out) and the fuzzy "AND" of the degree of memberships. 
A sum of these products is computed over all of the fired rules. The total is normalized by the sum of the minimum degree of 
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membership values for each rule fired. The output of the fuzzy processor represents the crisp estimate of the confidence in 
the velocity vector pair, ranging from 0 (low confidence) to 1 (high confidence). 

For example, given two pairs of vectors formed from two initial points competing for the same second particle as 
shown in figure 2 and in table 2: 



Figure 2: Two pairs of interacting vectors with separate initial points. Dashed line 
indicates incorrect vectors, solid line denotes correct vectors. 
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Table 2: Sample data for an interacting pair of 
vectors. 

We find the following combinations, and their respective input measures to the membership functions: 



Sep 

Mag 

MagDif 

Delta 

v,v. 

0 

7.3 

2.5 

25 

V,v 4 

4.3 

7.8 

1.4 

2.5 

V,v, 

1.7 

8.5 

5 

71.6 


2.6 

9.0 

3.9 

24.4 


Confidence 

0,24 

0.67 

0.11 

0.24 


Table 3: Fuzzy processor input measures and outputs for sample vector pair 
combinations 
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Looking at the first case, VjV 3 has the fuzzy set memberships: Sep (tts mall =0.875, 1X^-0, ; 

Mag {Ms^^.276, 1^=0.724, =0} ; MagDif (1^=0.504, ^=0, 1^=0}; Delta {^=0, ^=0, ^=1). 

With these fuzzy sets the following rules are fired: 

IF(Sep=Small AND Mag=Small AND MagDif=SmaU AND Delta=Large) THEN 
Conf_Out 1 =Med 

P_0Ut 1 =MIN(p SeptSrIull , M-Mig=Sm»ll’ ^MigDif-SnuIl’ l*Dehj=L»ige) 

END IF 

IF(Sep=Small AND Mag=Med AND MagDif=SmaU AND Delta=Large) THEN 
Conf_Out 2 =Low 

IL-OUtj— |X^jjg = Med’ M'MagDif^Small’ ^Delu=L»rge) 

END IF 

The fuzzy "AND" of the memberships in rule #1 gives p_outj=0.276 and rule #2 gives p_out2=0.504. The fuzzy processor 
output for this vector pair is given by: 


"m 

^ Conf _ Out i ■ |i_ out ( - 

Confidence = — (1) 

N i f 
i=l 


Confidence = 


0.5 0.276 + 0.1- 0.504 
0.276 + 0.504 


0.24 


( 2 ) 


where N R is the number of fired rules. The fuzzy processor output confidence level in each vector pair is shown in table 3. 
For the vector pairs listed, the V : V 4 pair has the highest confidence level of all the combinations. Initially, all of the vectors 
are assigned a confidence level of 0. The computed confidence level for the vector pair is compared against the current 
confidence level for each vector in the pair, and the maximum confidence value is stored for each vector. When all of the 
interacting vector pairs have been analyzed, the list of candidate vectors for each initial point will have confidence levels 
relative to one another. The vector in the list with the highest confidence level is assumed to be the most probable vector for 
the current initial point, and moved to the top of the list. 

A second pass through the data is performed to analyze non-interacting velocity vectors lists. An non-interacting 
vector is defined as a vector with multiple candidate vectors for a single initial point, but does not share a common second 
particle centroid with any other vector. The overall flow field pattern is used to determine the correctness of these outlying 
velocity vectors. The fuzzy processor is used to compare each vector in the list to all of the other previously analyzed velocity 
vectors in the flow field. Again, the vector with the highest confidence value in the list is assumed to be the correct vector. 

A final pass through the data is performed in which the highest confidence velocity vectors for each initial particle 
centroid are compared. All of the confidence levels are reset to 0. Each remaining valid vector is compared to all other valid 
vectors in the flow field. The resulting confidence level for each vector is then compared to a fixed threshold level. Only 
vectors with confidence levels above the set threshold level are considered valid. 


5. EXPERIMENT 

The experimental setup is shown in figure 3 and discussed in detail in reference 9. The test chamber consists of a 
10 cm diameter by 5 cm deep cylinder bored into a 12.7 cm square by 6.4 cm high plexiglass block. The reservoir was filled 
with 10 centistoke silicone oil and seeded with 12 pm pliolite particles. The flow was driven by the tubular heater in the 
center of the reservoir. The data were recorded onto a VCR using a gated intensified CCD camera. A large memory buffer 
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framegrabber was used to digitize the signals from the VCR at 640x480 pixels by 8 bits. Two successive video fields were 
acquired from a single camera for this low velocity flow. 



Figure 3: Schematic layout of experimental setup 
6. RESULTS AND DISCUSSION 

The flow field was allowed to reach steady state before the particle image data were acquired. The time interval 
between the two video fields was 1 second. The images were processed with a background threshold level of 10 grey levels. 
The particle centroids were determined via an intensity weighted mean calculation of the identified particle images in the 
field. The centroid processed images contained 316 and 318 particle images, respectively. The centroid error for this 
experimental configuration has been previously determined to be ±0.25 pixels. 10 


The information input into the fuzzy processor is a list of candidate velocity vectors for each initial particle centroid. 
The number of interacting velocity vectors is affected by the data density and the search region size R s . The 1 14 1 identified 
candidate velocity vectors are shown in figure 4. 

The output from the fuzzy processor is shown in figure 5. The resulting vector field map contains 209 velocity 
vectors. The total processing time was 7.6 seconds. The success rate of vector identification to number of particle image 
centroids in the image is 66 %. The minimum measured velocity vector displacement was 0.25 pixels and the maximum 
measured displacement was 17.7 pixels. The corresponding minimum and maximum velocities were 0.04 and 2.87 mm/s 
respectively. The measured velocity dynamic range was 70: 1 . There are 4 incorrectly identified velocity vectors shown in 
figure 5. The region near the fluid surface close to the heater is a region of high acceleration. The particle displacements in 
this region are on the order of or greater than 20 pixels. Many of the unsuccessful matchings have occurred in this region. 
The higher velocity region of the flow could be examined by decreasing the interframe acquisition time. 

Fuzzy PTV is similar to cross-correlation techniques used in normal PIV data reduction. Seeding densities 
comparable to those used for correlation data reduction yield good results. The local displacement information is used to 
determine the correct vector displacement. However, individual velocity vectors are detected, not the average displacement 
over an interrogation region. Hence, Fuzzy PTV can have better spatial resolution than correlation techniques. PTV 
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Future improvements to the fuzzy processor will be the inclusion of particle image shape and amplitude information 
as additional input measures. Inclusion of the particle image attributes further improves the robustness of the fuzzy tracking 
processor. The individual particle image shape and amplitude information can also be used to estimate the detected velocity 
vector accuracy as described in reference 10. 


7. CONCLUSIONS 

The application of fuzzy logic principles to the task of tracking particles recorded on two separate image frames has 
been demonstrated. The rule base developed is very simple and the fuzzy processor executed in software on a PC with good 
performance. Fuzzy techniques are simple to apply since no model of the flow field is required. Once the processor is 
established, no training is required as in neural network techniques. 

Fuzzy PTV is an attractive alternative to correlation techniques since no specialized processing hardware is required 
for data reduction or analysis. Direct digital recording simplifies the data processing. Fuzzy PTV requires only two 
independent image frames compared to traditional PTV techniques which require a sequence of images. Further 
enhancements to the fuzzy processor will increase its robustness and provide for individual velocity vector error estimates. 
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